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Abstract. I give a brief review of how X-rays from nearby galaxies are used 
as direct tracers of recent star formation. This leads to the conclusion that it 
is the most luminous point-like sources that are associated with star formation 
and that the majority of these are high-mass X-ray binaries. I then discuss a 
recent study that shows that ULXs are preferentially found in regions as young 
as or younger than typical Hii regions in their host galaxies. Finally, I describe 
a new study that attempts to determine the maximum luminosity of ULXs in 
the local universe by searching for them in interacting galaxies where the star 
formation rate is high. 



1 X-ray Sources as Tracers of Recent Star Formation 



David. Jones. &: Forman ( 19921 ) performed the first large survey of galaxies that 
showed a correlation between the total X-ray luminosity and galaxy-wide star- 
formation rate as measured by FIR (IRAS) luminosities. Although the X-ray 
data were not of a quality sufficient to distinguish among different types of X- 
ray sources, they argued on energetic grounds that most of the X-ray light came 
from high-mass X-ray binaries, with typical X-ray luminosities on the order of 
10^^ ergs s~^ and lifetimes of 10^ to 10^ yr rather than OB-star winds and SNe 
(bot h of which are mu c h fain t er though late stars liv e longer than XRBs; see 
also iHelfand &: Moraii (|200in : lDalton SarazinI (|1995I 1). 

The next major step forward awaited Chandra's superb angular resolution 
and broadband (CCD) spectral senstiv ity to isolate i ndividual X-ray sources in 
nearby galaxies. The definitive work is Grimm et al.l ([2003. ) . They showed that 
the point-source X-ray luminosity functions (XLFs) of individual normal and 
starburst galaxies all had roughly the same power law index over a wide range 
of source luminosities from <10^^ ergs s~^ up to about the highest luminosities 
observed of >10'^'' ergs s^^ (Figure 1, left panel). Furthermore, they showed that 
the overall normalization of these XLFs scaled with the total star-formation rates 
(SFRs) of the individual host galaxies. Thus, by rescaling observed XLFs to their 
universal XLF (Figure 1, right panel) one could deduce the star formation rate 
of the observed galaxy. Since the total X-ray luminosit y is just the integral of the 
XLF, it follows (rather circularly back to the work of iDavid. Jones. &: Forman 
(|1992| )1 that the SFR can also be deduced from a measure of the total X-ray 
luminosity even when individual sources cannot be resolved such as at high 
redshifts. 

Another key result of the iGrimm et al.l (120031 ^ study was that the (cumula- 
tive) XLF slope was rather flat, A^(> L) oc where A^(> L) is the number 
of sources with luminosity greater than L. This means that the most lumi- 
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Figure 1. Left: Cumulative XLFs for normal spiral and starburst galaxies. 
Right Same XLFs but normalized to the star-formation rate of the Antennae 
galaxies. No te that the sca l ed XL Fs occupy only a narrow range in iV(> L)—L 
space. From [Grimm "eTall (pool , by permission of Blackwell Publishing. 



2 Some Properties of Ultraluminous X-ray Sources 

The most X-ray luminous off-nucleus point-like objects in nearby galaxies are 
refered to as Ultraluminous X-ray sources or ULXs. Their high X-ray luminosi- 
ties (L ~^ 10'^^ ergs s~^ in the few tenths to ~10 keV range) imply high accretor 
masses; hence their interest as a possible evolutionary link between stellar-mass 
black holes (Mbh ~ 1 — 20 Mq) known throughout the Local Group and super- 
massive black holes (Mbh ^ 10^ — 10^ Mq) ubiquitous in the nuclei of galaxies. 

Numerous X-ray spectral and time- variability studies have shown that ULXs 
class, not steady thermal sources as expected from hot gas or SNRs but 
must be some form of X-ray binary (XRB). The association of the most lumi- 
nous sources with galaxy-wide star-formation indicators, such as FIR luminosity 
as discussed above, favor high-mass rather than low-mass XRBs. 

High-mass XRBs in our Galaxy and the Magellanic Clouds are those with 
short-lived, O or B type mass-donor stars. If this is the case also for ULXs, 
then they should be found preferentially in young star-forming regions such as 
Hll regions. Population synthesis models have shown that the highest X-ray 
luminosities, corresponding to the highest mass-transfer rates, are expected for 
high-mass XRBs aged roughly 4 to 10 Myrs old when the mass-donor compan- 
ion evolves off the main sequence. However, unambigously identifying optical 
counterparts to ULXs requires monitoring campaigns at high spatial resolution 
using Hubble or ground-based 8-m class telescopes. M ore fruitful are st udies of 
CMDs of objects in the neighborhood of ULXs, (e.g., ISoria et al.ll2005l l. which 



can narrow down the age of all likely counterparts in the region. 
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What is needed is a broadly-applicable method to easily quantify the star- 
forming properties of the local environments of ULXs without the expense of 
deep pointed obervations of individual ULXs. We have recently completed a 
survey using such a method ( Swartz et al.ll20()9l ). calibrated on the SDSS optical 
colors of known Hll regions, to distinguish young, star-forming regions in nearby 
galaxies. We have used these criteria to examine the 100x100 pc'^ regions around 
ULX candidates to determine if they are associated with the young stellar popu- 
lation. ULX environments in our study are designated as "star-forming" if they 
are as blue as, or bluer than a typical Hii region in its host galaxy (taking also 
into account the reddening effects of dust and of old underlying stellar popula- 
tions). A useful feature of this method is that colors are distance-independent 
quantities, and are measured on distance-independent regions. 

The 58 galaxies included in this study are a subset of a complete sample 
of 140 galaxies selected for X-ray analysis of their ULX populations. The X- 
ray sample galaxies are those that are in both the set of all galaxies within 
15 Mpc contained in the UGC with photographic magnitude nip < 14.5 mag 
and the set of all galaxies listed in the IRAS catalogs with a flux /fir > 
2Q-13.3 gj.gg QT^-'i g-i_ The 58 galaxies included in our study were those within 
the SDSS DR6 footprint with inclination i < 65°. 

We found that fully 60% of those ULX regions with sufficient signal (21/35) 
are in star-forming regions. By definition, these regions are as blue or bluer than 
typical H ii regions of their host galaxies and therefore are likely of an age also 
typical of Hii regions which is ^10 Myr; the characteristic lifetime of the least- 
massive LyC-producing stars, i.e., late O to early B stars of about 15 to 20 Mq. 

At the same time, we found that many of the most luminous ULXs in 
our sample are located in faint or non-starforming regions. Some of these may 
be heavily reddened by dust. However, we speculate that the most luminous 
ULXs (or equivalently, the most common phases of very high mass transfer) are 
biased towards early B-type donors with an initial mass of ~10— 15 Mq and an 
age ~10— 20 Myr (perhaps at the stage where the B star expands to become a 
blue supergiant). In that case, we expect very little residual Ha emission from 
their surroundings, the O stars having already evolved and died. Thus, we are 
beginning to constrain the age of the donor stars in ULX binary systems. We 
hope to soon constrain the age and mass of the accretor. 



3 A Search for the Most Luminous ULXs 

It was shown in S|T] that the X-ray luminosity function (XLF) of sources in 
the Milky Way and nearby normal spiral and starburst g alaxies is a smooth 
power-law distribution from -C 10^ ^ up to ~10^'' ergs s~^ ([Grimm et al. 20031 : 



Swartz et alll2004l : iLiu et al.ll2006l ). There are too few ULX candidates above 
~2xl0^'' ergs s~^ to state definitively whether or not the smooth distribution 
continues to higher luminosities. There are hints that a break or cutoff occurs 
at this point which, if real, imposes fundamental constraints on the physical 
properties of ULXs. 

The most important of these constraints is the mass of the accretor. We 
know empirically, from studies of Galactic sources, that the most reliable indirect 
mass indicator is the Eddington limit. If there is a cutoff in the luminosity 
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distribution at ~2xl0^'' ergs s~^, then this suggests a maximum black hole mass 
of ~100 Mq. The precise ma ss lirnit wo uld depend on the degree of anisotropy 
in the emission pattern (e.g., King 20081 ) and on the degree of super-Eddington 
luminosity; perhaps allowing for a factor of 1 to 2 lower limiting mass. On 
the other hand, if ULX luminosities extend beyond ~5xlO^'^ ergs (or if the 
luminosity function flattens), then masses >100 Mq are indicated, at least for 
the brightest sources. 

This is a critical dividing line. Theoretica lly, stellar cores up to ^70 M,^ can 
form equal mass black holes by direct collapse ( Yungelson et al.ll2008l : lHeger et al" 



20031 1 The final black hole mass may be even slightly larger if there is residual 



fall-back of any remaining stellar envelope. Above this core mass, stars disrupt 
via pair instability supernovae, leaving no remnant. Thus, something beyond 
normal single-star evolution will be required if ULXs are found to have lumi- 
nosities much higher than a few 10^" ergs s~^. 

Figure 2 shows the differen tial XLF for spiral galaxies taken from our orig- 
inal survey ( Swartz et al.|[2003 ). There are 57 spiral galaxies in this sample and 
97 ULX candidates but only a dozen above lO^'' ergs s~^. The fitted curves 
indicate a potential cutoff in the XLF at 2.07x10^'' ergs s~^ but the change in 
the fit statistic corresponds to an improvement over a single power law model 
at only 85% confidence. 




Figure 2. Differential XLF for spiral galaxies taken from ISwartz et al.l 
(2004). Three log-linear bins were used per decade in intrinsic 0.5-8.0 keV 
luminosity for the 97 ULX candidates. The best fit power law model and best 
fit exponentially cutoff power law model are shown. The exponentially cutoff 
power law is a better fit at 85% confidence. The upper curve shows the actual 
number of objects in each bin. 



We can estimate the number of ULXs needed to distinguish with high confi- 
dence between the power law and exponentially cutoff power law from the trends 
in the existing data. Assuming that the exponential cutoff power law is the true 



ULXs in Interacting Galaxies 



5 



functional form, we made a Monte Carlo simulation that computed the expected 
number of sources in each luminosity bin then generated a random variable with 
mean equal to the expected number and with a Poisson distribution. We then 
fit the simulated XLF with both models and compared the difference in the 
(maximum-likelihood) fit statistic. We conclude from this that about twice the 
current number of ULXs are needed to distinguish between these models at 
a >95% confidence level. This corresponds to ~30 additional ULXs above a 
luminosity 5x10^^ ergs s~^. 

Since both the number and peak luminosity of ULXs strongly correlates 
with galaxy-wide SFR and anti-correlates with nearest-neighbor distance, a log- 
ical choice of targets for our search are the high SFR interacting galaxies. Hap- 
pily, ISmith et al. ( 20071 ) have performed an in-depth study of the properties of 



a sample of pre-merger interacting galaxy pairs based on Spitzer luminosities 
and colors. They selected isolated, t idally disturbed, binary systems from the 
Atlas of Peculiar Galaxies ( ArpI 19661 ) that are within '^150 Mpc and that have 



large angular size to allow for good spatial resolution with Spitzer. Of the 35 
Arp pairs that met their criteria, we selected 7 galaxy pairs for a snapshot X-ray 
survey with Chandra which we believe will provide us with about 30 ULXs more 
luminous than 5x10^^ ergs s~^. These are being observed during the current 
Chandra observing round and the resul ts will be co r nbined with archival data 
from several other galaxy pairs in the Smith et al. ( 20071 ). So far, four Arp 



galaxy pairs have been observed and 13 ULX candidates have been detected 
above 5x10^^ ergs s~^. 
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